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We estimate the local number density of sources of ultra-high-energy cosmic rays 
^ ■ (UHECRs) based on the statistical features of their arrival direction distribution. 

CN , We calculate the arrival distributions of protons above 10^^ eV taking into account 

_!. \ their propagation process in the Galactic magnetic field and a structured intergalac- 

ff^ ■ tic magnetic field, and statistically compare those with the observational result of 

[■^^ . the Pierre Auger Observatory. The anisotropy in the arrival distribution at the 

^D \ highest energies enables us to estimate the number density of UHECR sources as 

10~ Mpc~ assuming the persistent activity of UHECR sources. We compare the 
estimated number density of UHECR sources with the number densities of known 
astrophysical objects. This estimated number density is consistent with the num- 
KN ' ber density of Fanaroff-Reily I galaxies. We also discuss the reproducability of the 

j^ ■ observed isotropy in the arrival distribution above 10 eV. We find that the esti- 

mated source model cannot reproduce the observed isotropy. However, the observed 
isotropy can be reproduced with the number density of 10~^-10~'^ Mpc~ . This fact 
indicates the existence of UHECR sources with a maximum acceleration energy of 
~ 10^^ eV whose number density is an order of magnitude more than that injecting 
the highest energy cosmic rays. 
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1 Introduction 



The origin of ultra-liigli-energy cosmic rays (UHECRs) has been highly un- 
known in spite of prolonged effort to construct larger UHECR observatories 
and to detect more events [1]. In cosmic ray spectrum, a sharply spectral 
steepening at around lO^^eV has been predicted theoretically by interactions 
with photopion production with cosmic microwave background (CMB) pho- 
tons, known as Greisen-Zatsepin-Kuz'min (GZK) steepening [2,3]. The feature 
of this steepening observed by the High Resolution Fly's Eye (HiRes) [4,5] and 
the Pierre Auger Observatory (PAO) [6,7] implies that astrophysical sources 
are much more dominant than top-down sources (for a review, see Ref. [8]) at 
the highest energies, though physical reasons for the extension of the energy 
spectrum beyond the GZK energies reported by the Akeno Giant Air Shower 
Array (AGASA) have been not understood yet [9,10]. Several objects to ac- 
celerate particles up to 10^'^ eV have been suggested, but there has been little 
observational evidence which is a UHECR source (see Ref. [8,11] and reference 
therein). 

The arrival distribution of UHECRs has information on the distribution of 
their sources, also including that on intergalactic and the Galactic magnetic 
field (IGMF and GMF). The AGASA and PAO reported statistically signif- 
icant anisotropy at small angular scale in the arrival distribution [12,13,14]. 
The small-scale anisotropy implies point-like sources, and enables us to con- 
strain the number of nearby UHECR sources. The source number density, n^, 
is one of important parameters to investigate the nature of UHECR sources. 
Comparing the estimated number density with the number densities of known 
astrophysical objects, we can constrain the object-type of UHECR sources. 
Several authors have estimated as Ug ~ 10~^-10~^ Mpc~^ using the published 
AGASA data above 4 x 10^^ eV assuming the same injection rate over all 
sources [15,16,17,18,19], and we also constrained the source number density 
as 10"^ Mpc~^, assuming a model that the injection rate is proportional to 
the luminosity of galaxies [19]. 

The PAO reported the correlation between the arrival directions of UHECR 
above 5.7 x 10^^ eV and the positions of nearby active galactic nuclei (AGNs) 
listed in the 12th Veron-Cetty & Veron catalog [20] within the angular scale of 
3.1° [13,14]. However, most of the PAO-correlated AGNs are Seyfert galaxies 
and LINERs, which have much weaker activity than radio-loud AGNs like 
Fanaroff-Reily H (FR H) galaxies [21]. To understand whether these galaxies 
with weak activity are really UHECR sources or not, another information on 
UHECR sources is required. 

The PAO estimated a lower limit of the UHECR source number as 61 based on 
anisotropy in the arrival distribution of their detected events, simply assuming 



Poisson statistics and not taking into account UHECR propagation [14] . This 
is certainly a lower limit of the number of the sources, but is not quantity which 
can be compared with astronomical observables because we do not know what 
radius they are included in. Their number density is an observable. Taking 
into account UHECR propagation, we can also estimate plausible value of it, 
not limit. 

In this study, we simulate the arrival distribution of UHECRs above 10^^ 
eV, taking into account their propagation process in intergalactic and the 
Galactic space. We extract an anisotropy signal from the simulated arrival 
distribution and then compare this signal with that observed by the PAO 
above 5.7 x 10^^ eV to estimate the source number density. The reports for 
the correlation by the PAO with nearby large-scale structure [13,14,22,23] 
show that it plays a crucial role on the arrival distribution of UHECRs. Thus, 
we adopt the models of the IGMF and source distribution which reproduce 
the local structure actually observed around the Milky Way, developed in our 
previous work [18]. We also discuss the isotropy in the arrival distribution at 
~ 10^^ eV and estimate the source number density for the lower energies. The 
composition is assumed to be pure protons. 

This paper is structured as follows. In section 2, our calculation method for 
calculating UHECR arrival distribution and a statistical method to estimate 
the source number density are explained. In section 3, the calculation results 
are shown and we discuss the results and conclude in section 4. 



2 Methods 



We estimate the number density of UHECR sources as follows. First, we cal- 
culate arrival distributions of UHE protons based on our source models with 
several number densities. Their propagation process in intergalactic and the 
Galactic space is taken into account. The number of simulated protons and 
the threshold of their energies are set to the same as those detected by the 
PAO. Next, auto-correlation functions, which are an indicator of small-scale 
anisotropy in the arrival distribution, are calculated from our simulated events. 
Finally, comparing these auto-correlation functions with that calculated from 
the observed events, we find the best-fit value of the source number density. 

We calculate the arrival distribution of UHE protons by a method used in Ref. 
[19]. A characteristic of this method is to adopt a structured IGMF model 
which reproduces the observed structure in local universe, developed in Ref. 
[18]. The strength of the IGMF is normalized at the center of Virgo cluster as 
B = 0.0, 0.1, 0.4 and l.OfiG. This method also enables us to take account of the 
deflection by GMF. A GMF model with bisymmetric spiral structure proposed 



by Ref. [24] and the same model parameters as Ref. [25] is adopted. The 
maximum distance of UHECR sources is set to IGpc, which is comparable with 
the energy-loss length of Bethe-Heitler pair creation with the CMB. This is 
sufficient to consider cosmic rays down to 10^^ eV. There are two improvements 
for this study. One is the angular resolution of UHECR experiments. In Ref. 
[19], since we consider the AGASA data, we adopt the angular resolution of 
~ 2°. In this study, 1° is adopted for the PAO events, which corresponds to the 
angular resolution of the PAO [14]. The other is to take the non-uniformity of 
the exposure of ground array into account. 

The exposure of ground array is not uniform because of the daily rotation 
of the earth. Since the variation in right ascension in a day can be neglected 
[14], the dependence of the exposure is simply written as a function of the 
declination of arriving cosmic rays, 6, [26] , 

uj{6) oc cos(ao) cos(5) sin(Q;m) + a^n sin(ao) sin(5) (1) 

where a^ is given by 

if e > 1 

ttm = < vr if ^ < -1 (2) 

cos~^(.^) otherwise. 



and 



^ cos{9m) - sin(ao) sin(^) 
cos(ao) cos((5) 



Here, oq is detector's terrestrial latitude, and 6„i is the maximum zenith angle 
for the experimental cut. For the PAO, Cq = 35.2° and 6 = 60° are adopted 

[26]. 

In this study, we adopt two source models. Both models are constructed from 
Infrared Astronomical Satellite Point Source Redshift Survey {IRAS PSCz) 
catalog of galaxies [27]. The IRAS catalog consists of 14,677 galaxies above 
0.6 Jy with redshift and it sky coverage is 84%. So it is the best catalog to 
reflect the local large-scale structure. We correct the selection effect and galax- 
ies in 16% of unobserved area using the luminosity function estimated by Ref. 
[28] , and then regard subsets of the corrected IRAS galaxies as UHECR source 
distributions. Each galaxy is adopted as a source with the probability propor- 
tional to its luminosity. We perform the source selection 100 times for every 
Us- One of the models is a source model with the same cosmic ray injection 
rate over all sources (called model A below), and the other is with the injection 



rate proportional to infrared luminosities of galaxies (called model B below). 
The maximum acceleration energy of protons is not dependent on luminosities 
of the sources. This corrected catalog is also used for the construction of our 
IGMF model. 

Auto-correlation function is one of good indicators of small-scale anisotropy 
of UHECR arrival distribution. Auto-correlation function is defined as 



where is the separation angle of a pair of events. The interval A^ is set to 
be 1°. In order to investigate the goodness of fit between the auto-correlation 
functions calculated from simulated events and observed one, we also define 
Xe^ax as 



X.6-n 



Nun\ 



n=0 Cr{(^n) 



where w{6) is an average of the auto-correlation function calculated from a 
source distribution and a{9) is the standard deviation of w{9) due to the 
finite number of simulated events. Random event selection is performed 100 
times for every source distribution. Wohs{0) is the auto-correlation function 
calculated from observed data, ^max is the maximum angular scale to consider 
the small-scale anisotropy, A^bin = ^max/A^ is the number of angular bins. 



3 Results 



3.1 Estimation of Us from the small-scale anisotropy 



We estimate UHECR source number density using an anisotropy signal in the 
27 events published by the PAO whose energies are above 5.7 x 10^^ eV [14]. 
First, we check the anisotropy signal. Fig. 1 shows auto-correlation function 
calculated from the 27 events {histogram) . The auto-correlation function pre- 
dicted from random distribution weighted by the exposure of the PAO is also 
shown. The error bars are the standard deviations due to the finite number of 
events. The event realization is performed 1000 times. We can see statistically 
significant anisotropy at small angular scale against isotropic distribution. We 
set 6'max = 5° in Eq. 5 to extract the small-scale anisotropy. If ^max ~ 10° is 
adopted, main results below are unchanged. 



Next, we calculate auto-correlation functions from our simulated arrival dis- 
tributions of UHE protons, which are compared with the histogram in Fig. 1 
to investigate the number density of the sources which can best reproduce the 
small-scale anisotropy observed by the PAO. If Ug is extremely small, strong 
anisotropy is expected because only few sources can contribute to the arriv- 
ing cosmic rays. On the other hand, in the case of much many sources, the 
arrival distribution is expected to be close to isotropy. Thus, there should be 
the best-fit value of the number density. 

Fig. 2 shows X5 calculated for the different number densities of UHECR 
sources. The normalization of the IGMF strength considered is 0.0 {crosses), 
0.1 {triangles), 0.4 {squares), and 1.0 /xG {pentagons) respectively. The error 
bars represents the standard deviation of Xs? which is estimated from 100 re- 
alizations of the source distribution. The plots at the same number density 
are a little shifted horizontally for visibility. Note that UHECR sources within 
5 Mpc are artificially neglected in this calculation to reproduce isotropy ob- 
served at lower energies (~ 10^^ eV). This point will be discussed in section 
3.2. Such nearby sources also predict anisotropy stronger than observed on at 
the highest energies. 

In the left figure, xs calculated in the model A is shown. The GMF is consid- 
ered in the lower panel while not in the upper panel. For n^ ~ 10^^ Mpc~^, 
there is zero or one source within 100 Mpc. This extremely few number of 
sources leads to anisotropy much stronger than the observational result. Thus, 
rig ~ 10^^ Mpc^^ is too small to reproduce the observed arrival distribution. If, 
on the other hand, n^ ~ 10^^ Mpc^^, which is comparable with the number of 
bright galaxies, high-precision isotropy of the arrival distribution is realized. 
This another extreme case cannot also reproduce the observed anisotropy. 
Thus, X5 ^^ ci function of Ug has the minimum at the intermediate number 
density. In the upper panel, Xb is minimized at around Ug ~ 10~^-10~^ Mpc^^, 
almost independent of the IGMF strength in our IGMF model. In our IGMF 
model, since about 95% of volume within 100 Mpc has no magnetic field, the 
highest energy protons are defiected only in the neighbourhood of their sources 
where the universe is magnetized. If the universe is strongly magnetized, for 
example, uniformly with ~ lOOnG, the source number density is allowed to be 
smaller because strong defiections makes the arrival distribution to be isotropic 
even at the highest energies. In the lower panel, xs is minimized at the same 
number densities. The constrained number density is not dependent on the 
existence of the GMF since the clustering signal is almost not varied because 
of the coherence of the GMF. Note that the GMF can changes the arrival 
directions of UHECRs efficiently and positional correlation with their sources 
as shown in Ref. [29]. 

We also estimate the source number density in the model B by similar discus- 
sion above. In the right figure, Xs iii the model B is shown, xs is minimized 



at around Us ~ 10~^-10~^ Mpc~^, which is about an order of magnitude more 
than in the model A. In the model B, luminous sources strongly contribute to 
the arriving cosmic rays and, on the other hand, weaker sources do not almost 
contribute the flux of cosmic rays though these are counted as sources. Thus, 
the constrained number density is effectively larger. 

The number densities estimated from the anisotropy signal of the PAO based 
on the two source models (A & B) are consistent with the results of our 
previous study using the AGASA data [19]. The uncertainty of about one 
order of magnitude will be reduced by 5 years observation by the PAO, as 
pointed out in Ref. [19]. 

The model A and B are two extreme cases. In the model A, all sources are 
identical. This model can be well applied to astrophysical objects which have 
hardly personality regarding emitted energy, or strongly active sources like 
FR II galaxies and Gamma Ray Bursts (GRBs) (though GRBs are transient 
sources). On the other hand, the model B is applied to objects common in 
the universe like bright galaxies. When the model B was constructed, the 
IRAS galaxies with luminosities of IO'^-IO^^Lq were adopted, where Lq is 
the Solar luminosity, 3.826 x 10^^ W. The luminosities are over 5 orders of 
magnitude. The width of magnitude is realized by less active objects like 
bright galaxies and Seyfert galaxies. In fact, UHECR sources have intermediate 
nature between the two models. Assumed to be highly active objects, UHECR 
sources has the number density close to that in the model A. Thus, the number 
density of UHECR sources can be estimated as ~ 10^^ Mpc~^. 

Compared with the number densities of several candidates of UHECR sources, 
this estimated value can constrain the object-type of UHECR sources. Table 
1 shows local number densities of several astrophysical objects. The objects 
which have smaller number density than the estimated one are not main con- 
tributors of the arrival highest energy cosmic rays. Both FR II galaxies and BL 
Lac objects have been plausible UHECR sources theoretically [30,11], but they 
are not mainly contribute to the flux of UHECRs. Note that this constraint 
does not reject them as UHECR sources. They are not main contributors. The 
number density of bright galaxies (deflned as —22 < Mbj. < —18, where Mbj. 
is absolute magnitude in B-band [31]) is much larger than our constraint. The 
generation of the highest energy cosmic rays is not common phenomena in the 
local universe. The number density of Seyfert galaxies is an order of magni- 
tude larger than the constrained value, and therefore it is allowed if ~ 10% 
of such Seyfert galaxies has activity to accelerate UHECRs. However, Seyfert 
galaxies are generally not expected to be UHECR sources. For example, the 
PAO-correlated AGNs, most of all is Seyfert or LINER, do not show signif- 
icant jet activity [21]. Therefore, there are no reasons for expecting them to 
accelerate cosmic rays up to the highest energies at all in the jet paradigm. 
The number density of FR I galaxies is ~ 10^^ Mpc"'^ [34], consistent with the 



constraint. FR I galaxy is one of plausible acceleration site up to the highest 
energies. Centaurus A (Cen A), the nearest PAO-correlated AGNs, is classi- 
fied into FR I galaxy. At a hot spot in the jet in this galaxy, an estimation of 
the maximum acceleration energy of protons is ~ 10^'^'^ eV [11]. If the other 
FR I galaxies have similar configurations, it is possible that FR I galaxies are 
main contributors of the highest energy cosmic rays. 



3.2 Isotropy at around 10 eV 



At the previous section, we artificially neglect sources within 5 Mpc from the 
Galaxy. This is because very nearby sources inevitably produce anisotropy in 
arrival distribution of UHECRs at around lO^^eV. Since all observatories on 
UHECRs have reported isotropy in that energy range [12,35,36]. it is inconsis- 
tent with the observational results. We discuss the isotropy in this subsection. 

Fig. 3 shows the arrival distributions of protons above 10^^ eV simulated from 
a source distribution in the model A with Ug ~ 10~^ Mpc~^, which is about 
one order of magnitude larger than the number density constrained in the 
previous section. In the left figure, the sources within 5 Mpc are neglected. 
The number of events is set to 1672, which corresponds to that detected by 
the PAO [36]. The IGMF strength is fi = 0.4^G. 

In the right figure, strong anisotropies appear at the positions of (£, b) ~ 
(—50°, 20°) and (100°,— 85°), which are generated by protons injected from 
Cen A and NGC 253, a famous starburst galaxy, respectively. These distances 
are 4.6 Mpc and 4.2 Mpc in the IRAS PSCz catalog, respectively [27]. Note 
that the distance of Cen A is different among many catalogs. If these nearby 
objects are really UHECR sources, strong anisotropy is inevitably predicted. 
The PAO reported that the arrival distribution of UHECRs above 10^^ eV is 
consistent with isotropic distribution with 95% confidence level [36]. Thus, the 
anisotropy signals are inconsistent with the PAO results. Such anisotropy is 
predicted as long as there are UHECR sources in the nearby universe even if 
Hg ~ 10~^ Mpc~^ is adopted. On the other hand, the anisotropies disappears 
in the left figure. This fact is also true in the model B, as shown in Fig. 4. 
Thus, these nearby objects are not UHECR sources. 

In fact, the observed high-level isotropy cannot be easily reproduced by ne- 
glecting sources within 5 Mpc. In order to see this, we compare the auto- 
correlation function of simulated events above 10^^ eV with that predicted 
from isotropic distribution. Table 2 shows the fraction of the number of source 
distributions which can reproduce the arrival distribution consistent with the 
isotropy within some errors written in the table, a means the standard de- 
viation of the auto-correlation function calculated from random distribution. 



which is due to the finite number of events. The normahzation of the IGMF 
strength considered is 0.1 and 0.4 /xG. In the model A, there are few source 
distributions with Ug ~ 10^^ Mpc~^ which can reproduce the isotropy within 
2a. Even for Ug ~ 10~^ Mpc~^, it is less than 10 though isotropy is found by 
eye in Figs. 3 and 4. If the number density is ~ 10^^ Mpc~^, the isotropy can 
be well reproduced. 

In addition to the structured IGMF, we take account of a uniform turbu- 
lent magnetic field with the strength and coherent length of -Bturj h, re- 
spectively. Faraday rotation measurements of distant quasars shows BlJ"'"^ < 
(1 nG)(l Mpc)^/^ [37]. According to this constraint, we adopt i?tur = 1 nG 
and /c = 1 Mpc. The results are also shown in Table 2. We are able to find 
that the number of the source distributions which can reproduce the isotropy 
is increased, compared with that without the turbulent field. However, for the 
estimated number density, the number is still too small to naturally reproduce 
the isotropy. When the source number density an order of magnitude larger 
than the estimated one, the number is dramatically increased to 70-80% within 
3cr in the model A. If the number density is ~ 10~^ Mpc~^, the isotropy can 
be well reproduced. 

It is much larger than the number density estimated based on the small- 
scale anisotropy at the highest energies. This fact indicates the existence of 
sources which contributes to lower energy cosmic rays, whose number density 
is 10~^-10~^ Mpc~^. In other words, there are many sources with the maximum 
acceleration energy of ~ 10^^ eV. Note that the uniform IGMF does almost 
not affects the estimation of the number density at the highest energy. 



4 Discussion & Conclusion 



In this study, we estimated the number density of UHECR sources based on 
the statistical features of the arrival distribution of UHECRs observed by the 
PAO. We simulated the arrival distributions of protons above 10^^ eV, taking 
into consideration their propagation process in the Galactic and intergalactic 
space. The IGMF model adopted was associated with the observed large- 
scale structure. Comparing the simulated arrival distributions with the data 
observed by the PAO statistically, we estimated the number density which can 
best reproduce the observational result. The anisotropy signal above 5.7 x 10^^ 
eV led to n^ ~ 10~^ Mpc~^ which is consistent with that of FR I galaxies, 
which are a candidate to accelerate protons up to 10^" eV. We also focused 
on isotropy in the arrival distribution at around 10^^ eV, and then found 
Us ~ 10~^-10~^ Mpc~^, which is one or two order of magnitude larger than 
that estimated from the anisotropy. 



In this calculation, cosmological evolution of UHECR sources is not consid- 
ered. Since protons with ~ 10^^ eV can reach the earth from sources at the 
distance of IGpc {z ~ 0.25), which is comparable with the energy-loss length 
of such protons by Bethe-Heitler pair creation, it might be possible that the 
cosmological evolution could reproduce the isotropy even if the local number 
density of UHECR sources is comparable with ~ 10"'^ Mpc^^. However, a typ- 
ical evolution factor is (1 + z)^, and therefore the number density at z ~ 0.25 
is only twice more than the local one. Such a small factor could not change 
the estimated number density by two order of magnitude less. The difference 
between the two estimated number densities is significant. 

The difference can be interpreted as the evidence of the existence of UHECR 
sources which can accelerate protons up to ~ 10^^ eV assuming the persistent 
activity of UHECR sources, n^ ~ 10~^-10~^ Mpc~^ is comparable with the 
number densities of bright galaxies or Seyfert galaxies. The proton accelera- 
tion up to 10^^ eV might be common in the universe. We could also interpret 
this difference as transient activity to emit UHECRs. If bursting sources, like 
GRBs, are assumed, apparent number density, which corresponds to Ug esti- 
mated in this study, depends on the threshold of UHECR energies because the 
dispersion of the time-delay is larger at lower energies [38]. However, quanti- 
tative discussion on this possibility exceeds the scope of this study. It is next 
study of ours. 

In order to reproduce isotropy in the arrival distribution of UHECRs above 
10^^ eV, sources within 5 Mpc were artificially neglected. In these sources, Cen 
A, the nearest radio-loud AGN, is involved. The isotropy at around 10^^ eV 
implies that it is not UHECR sources in our persistent source model. However, 
radio-loud AGNs have been plausible site for particle acceleration up to the 
highest energy [11], and the two of the highest energy events of the PAO are 
correlated with the position of Cen A within ~ 3°. Whether Cen A is really a 
UHECR source or not is a key to understand the mechanism to generate the 
highest energy cosmic rays. 

If Cen A is not a UHECR source, we can interpret that the source of cosmic 
rays arriving in the direction of Centaurus is behind Cen A and more distant, 
so that the strong anisotropy is not generated. The 2 events towards Cen A 
are also positionally correlated with NGC 5090 with ~ 3°, a radio galaxy with 
the distance of ~ 40 Mpc. This might be a real source though the activity 
is weak like the other PAO-correlated AGNs. These events are also towards 
Centaurus cluster whose distance is about 40 Mpc. An idea to generate UHE 
particles is cluster accretion shock [39,40,41]. However, this scenario can ac- 
celerate protons up to ~ 10^^ eV. 

If Cen A is a UHECR source, we propose several possibilities not to gener- 
ate the anisotropy. One is UHECR composition. In the case of heavy dom- 
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inated composition at energies above 10^^ eV, the trajectories of UHECRs 
are deflected Z times more than those of protons, and the discussion on the 
isotropy in this paper is not apphed. Several composition measurements imply 
the existence of some fraction of heavy elements for all large uncertainty on 
hadronic interaction models in extensive air shower [42,43]. The other is that 
the UHECR production is transient. For a bursting source, the energy of cos- 
mic rays observed at present is in narrow energy range since the time-delay of 
cosmic rays depends on their energies [38]. In this brief picture, the 2 events 
are arrived from Cen A, and the lower energies will come in the future. The 
detailed discussion on the two possibilities is our near future plan. 

The estimated number densities have some uncertainty because of the small 
number of detected events. The uncertainty to determine the number density 
at the highest energies can be reduced by increasing observed events [19] and 
significant anisotropy would be observed in the arrival distribution at ~ 10^^ 
eV in the near future [44,22]. The dramatically increase of detected events in 
the near future will provide us more useful information on UHECR sources. 
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Fig. 1. Auto-correlation functions calculated from the 27 PAO events (histograTn) 
and random distribution weighted by the exposure of the PAO. The error bars are 
the standard deviations due to the finite number of events. The event realization is 
performed 1000 times. 
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Fig. 2. Xb for different number densities of UHECR sources in the model A (Ze/t) 
and model B {right). The error bars originate from 100 times source selection. The 
GMF is considered in the lower panel while not in the upper panel. 
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Fig. 3. The arrival distributions of 1672 protons above 10^^ eV (= 10 EeV) calculated 
from a source distribution in the model A with n^ ~ 10~^ Mpc~^ in which the 
sources within 5 Mpc are included (right) and neglected (left). The exposure of the 
PAO is taken into account. The IGMF strength is B = 0.4/xG and the GMF is 
neglected. 
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Fig. 4. The same as Fig. 3, but in the model B. 



Table 1 

Local number densities of several active objects. 



Object 



Density [Mpc ] Reference 



Bright galaxy 
Seyfert galaxy 
Bright quasar 
Fanaroff-Reily 1 
Fanaroff-Reily 2 
BL Lac objects 



1.3 X 10^2 
1.25 X 10-3 

1.4 X 10"^ 
8 X 10-5 
3 X 10-*^ 
3 X lO-'^ 



Ref. [31] 
Ref. [31] 
Ref. [32] 
Ref. [34] 
Ref. [33] 
Ref. [33] 
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Table 2 

The numbers of the source distributions which satisfy isotropy in the observed 
arrival distribution of UHECRs above 10^^ eV within the error bars shown in table, 
in the case of the model A (without parentheses) and B (in parentheses). The total 
number of source distributions in each parameter set is 100. 
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